Distribution of Ni Yields of Type la Supernovae and 
Implications for Progenitors 

Bo Wang^'^*, Xiang-Cun Meng^'^, Xiao-Feng Wang^''^ and Zhan-Wen Han^ 
^ National Astronomical Observatories/ Yunnan Observatory, Chinese Academy of Sciences, 
Kunming 650011 

^ Graduate School of Chinese Academy of Sciences, Beijing 100049 
^ Tsinghua Center for Astrophysics (THCA) and Department of Physics, Tsinghua 
University, Beijing 100084 

Astronomy Department, University of California at Berkeley, CA 94720, USA 
22 June 2007 

Abstract The amount of ^^Ni produced in type la supernova (SN la) explosion 
is probably the most important physical parameter underlying the observed corre- 
lation of SN la luminosities with their light curves. Based on an empirical relation 
between the ^^Ni mass and the light curve parameter Ami5, we obtained rough 
estimates of the ^^Ni mass for a large sample of nearby SNe la with the aim of ex- 
ploring the diversity in SN la. We found that the derived ^^Ni masses for different 
SNe la could vary by a factor of ten (e.g., Mni — 0.1 — 1.3 Mq), which cannot be 
explained in terms of the standard Chandrasekhar-mass model (with a ^^Ni mass 
production of 0.4 - 0.8 Mq). Different explosion and/or progenitor models are 
clearly required for various SNe la, in particular, for those extremely nickel-poor 
and nickel-rich producers. The nickel-rich (with Mm > 0.8 Mq) SNe la are very 
luminous and may have massive progenitors exceeding the Chandrasekhar-mass 
limit since extra progenitor fuel is required to produce more ^^Ni to power the light 
curve. This is also consistent with the finding that the intrinsically bright SNe la 
prefer to occur in stellar environments of young and massive stars. For example, 
75% SNe la in spirals have Amis < 1-2 while this ratio is only 18% in E/SO galax- 
ies. On the other hand, the nickel-poor SNe la (with Mni < 0.2 Mq) may invoke 
the sub-Chandrasekhar model, as most of them were found in early-type E/SO 
galaxies dominated by the older and low-mass stellar populations. This indicates 
that SNe la in spiral and E/SO galaxies have progenitors of different properties. 
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1 INTRODUCTION 

It is well known that Type la Supernovae (SNe la) are excellent cosmological distance indicators 
due to their high luminosity and remarkable uniformity. The results of the High-Z Supernova 
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Search Team (Riess et al. 1998) and the Supernova Cosmology Project (Perlmutter et al. 1999) 
astonishingly showed that the expansion of the universe is accelerating. Although SNe la are 
very important objects in modern cosmology, several key issues related to the nature of their 
progenitor systems (Hachisu et al. 1996, 1999; Li & van den Heuvel 1997; Han & Podsiadlowski 
2004; Meng et al. 2006) and the physics of the explosion mechanisms (Hillebrandt & Niemeyer 
2000; Ropke & Hillebrandt 2005a) are still not understood. 

It is widely accepted that SNe la are thermonuclear explosions of carbon-oxygen (C-0) 
white dwarfs (WDs) accreting matter from their companions (for a review see Nomoto et al. 
1997). The energy released from the burning till nuclear statistical equilibrium (NSE) completely 
destroys the C-0 WD. The optical/IR light curves are powered by the radioactive decay of ^^Ni 
^ 56co ^ 56pe (Colgate & Mckee 1969). 

Over the last decades, it has been found that there exists spectroscopic diversity among 
SNe la. In an attempt to quantify the rate of spectroscopically peculiar SNe la in the existing 
observed sample, Branch et al. (1993) compiled a set of 84 SNe la and foimd that about 
83% — 89% of the sample are normal (or sometimes "Branch normals"). According to Li et al. 
(2001), however, only 64% ± 12% of observed SNe la are normal in a volume-limited search. 
The total rate of pecuhar SNe la could be as high as 36% ± 9%. The rates are 16% ± 7% 
and 20% ± 7% for SN 1991bg-nke events and SN 1991T-Hke objects, respectively. Moreover, 
some of the peculiar SNe la apparently deviate from the relation between the light-curve shape 
parameter Ami5 and the luminosity (e.g., Wang et al. 2006), which probably suggests that 
SN la luminosity is not a single parameter of the light curve shape. Note that Wang et al. 
(2005) presented a single post-maximum color parameter AC12 {B — V color ~ 12 days after 
the B-band light maximum), which empirically describes almost the full range of the observed 
SN la luminosities and gives tighter correlations with SN la luminosities, but the underpinning 
physics is still not understood. Stritzinger et al. (2006a) recently noticed that no explosion 
model has been proposed yet to successfully account for the factor of ten spanned by the range 
in the observed bolometric luminosity. 

In this paper, we aim to understand the origin of the progenitors and explosion models 
with a large well-observed sample of nearby SNe la. According to a simple empirical relation 
between the ^^Ni mass (Mn;) and the decline rate parameter Amis, we can estimate the mass 
of radioactive material ^^Ni synthesized in SN la explosions. By studying the distribution of 
such a parameter for SNe la, we hope to obtain some clues to the SN la diversities. Moreover, 
we also investigate the correlation of SNe la distribution with their host galaxies, which might 
also set constraints on the SN la progenitors. 

The paper is organized as follows. In Section 2, we provide a brief description of the ob- 
servational data examined in this study. Section 3 shows the distribution of the decline rate 
parameter Ami5, the relative radial distance rsN/-R25 and the derived ^^Ni mass. In particular, 
a simple empirical method is introduced and used to estimate the ^^Ni mass for a large number 
of well-observed SNe la. A discussion and conclusions are then given in Section 4. 

2 DATA 

In this paper, we consider a sample of 111 well-observed nearby SNe la with Z < 0.1, of these 
109 were compiled by Wang et al. (2006) (see tables 1 and 2 in their work), and SNe 2002cs and 
2002dj are from unpublished KAIT SN database. The sample (excepting SNe 1937C, 1972E 
and 1974G) includes only SNe with CCD measurements and those observed not later than 8 
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days after the -B-band light maximum and with more than five photometric points. A large 
fraction of these observations were obtained by the earUcr Calan/Tololo SN survey, the CfA 
I/II SN monitoring campaign, and the Las Campanas/CTIO observing campaign. 

The sample contains two important parameters (i.e.. Amis and rsN/.R25)- The decline rate 
parameter Amis was defined by Phillips (1993) as the amount in magnitudes that the S-band 
light curve decays in the first 15 days after maximum light, and the decline rate was corrected 
for a small reddening effect (Phillips et al. 1999). rsN/-R25 is the de-projected galactocentric 
distance of the SNe in their respective host galaxies in units of the galaxy radius R25, where rsN 
is the de-projected radial distance of the SN la from the galactic center and R25 is the apparent 
semimajor axes of the SN host galaxy isophote having /i^ = 25.0 mag arcscc"^. The 111 wcll- 
observcd sample contains 90 spectroscopically normal SNe la, 13 SN 1991T-like objects, six 
SN 1991bg-like events and two truly pecuUar SNe la, SNe 2001ay and 2002cx. The 111 sample 
can be divided into four groups by the morphological type of their host galaxies: 77 are in 
spiral galaxies, 28 in E/SO galaxies, two (SNe 1972E and 2000fa) in irregular galaxies, and four 
(SNe 1992bl, 1999aw, 1999ej and 2002cx) in galaxies of uncertain or intermediate types. Again, 
the sample of 111 contains 98 SNe la with relative radial distances while the remainder are 
uncertain. Out of the 98 SNe la with known distances, 26 are from E/SO galaxies, 70 from 
spiral galaxies, one (SN 1992bl) from a SO/a galaxy and one (SN 2000fa) from an irregular 
galaxy. 

3 ANALYSIS AND RESULTS 

In this section we investigate the distributions of Amis and the derived ^^Ni mass versus 
^sn/-R25. In particular, we provide a short overview of the method to estimate the ^^Ni masses 
for a large number of well-observed SNe la. 

3.1 Distribution of Ami5 

Using optical light curves of SNe la, Phillips (1993), Hamuy et al. (1996a) and Phillips et al. 
(1999) established a relationship between the absolute magnitude at S-band light maximum and 
the initial decline rate of the light curve. The characteristic parameter is Amis. When Amis 
is properly corrected for, the SNe la prove to be a high precision distance indicator, yielding 
relative distances with uncertainties ~ 7% - 10% (Hamuy et al. 1996b). Kasen & Woosley (2007) 
confirmed that the brighter SNe la have broader B-band light curves with smaller decline rates 
around optical maxima than the dimmer SNe la. Moreover, the determination of the decline rate 
parameter Amis is available for a large number of objects. By analyzing Amis, it is possible 
for us to peer into some intrinsic properties of SNe la, because it is basically independent of 
distance and reddening. 

Hamuy et al. (1995, 1996a) first noticed that brighter (also smaller Amis) SNe la are 
preferentially located in late-type galaxies for the nearby sample of SN la (see also Gallagher et 
al. 2005). In Figure 1, we show the distribution of Amis for spiral and E/SO galaxies, where the 
mean error for Amis is 0.07. It is obvious that the distribution of Amis is different in spiral and 
E/SO galaxies. A K-S test shows that the probability of these two distributions coming from the 
same population is as low as P = 3.3 x 10~^ (see also Altavilla et al. 2004; Delia Valle et al. 2005 
for similar argument). Although both SNe la in spiral and E/SO galaxies show a wide range of 
Amis, the fraction with Amis < 1-2 is 75% in spirals and only 18% in E/SO galaxies. The fact 



4 



B. Wang et al. 




Fig. 1 Distribution of Amis for 28 SNe 
la in E/SO galaxies (solid line) and 77 in 
spiral galaxies (dotted line). 



Fig. 2 Distribution of Amis for differ- 
ent SN la types. The dotted line shows 
the distribution of 90 normal SNc la. 
The cross-hatched area and the simply 
shaded area are for 13 SN 1991T-like events 
and 6 SN 1991bg-like objects, respectively. 
Moreover, the gray-colored area and the 
black area show SNe 200 lay and 2002cx, 
respectively. 



that faster decliners (i.e., fainter SNe la) dominate in E/SO galaxies is apparently related to the 

age of the progenitors since E/SO galaxies arc believed to have older low-mass stars. Recently, 
Sullivan et al. (2006) demonstrated that the stretch factor of SNe la (another form of Amis) in 
low- and high-redshift spiral galaxies bears the same distribution at a confidence level of 95%, 
while the distributions in elliptical galaxies are different as the faster decliners (also the fainter 
SNe la) are usually seldom observed at high redshift probably due to the selection effect and/or 
the SN evolution effect. Since the decline rate Amis is tightly correlated with the stretch factor 
s, this difference found in elliptical galaxies holds true for Amis. 

Figure 2 shows the distribution of Amis for different types of SNe la. Compared with 
normal SNe la, SN 1991T-like events have a flatter dechne rate of 0.95±0.10 while SN 1991bg- 
like objects have a much faster decline rate of 1.90±0.10. In the figure SNe 2001ay and 2002cx 
are truly peculiar events. SN 2001ay has the slowest decline rate in our sample, and Howell 
& Nugent (2004) suggested that SN 2001ay is a SN la with the broadest light curve. Li et al. 
(2003) noted that SN 2002cx has a pre-maximum spectrum like SN 1991T (the classical hot SN 
la), a luminosity like SN 1991bg (the classic subluminous event), a slow late-time decline, and 
unidentified spectral lines (which is probably not a SN la as it showed a spectroscopic feature 
extremely similar to those in some faint type II-P SNe at nebular phase, private communication 
with Dr. Weidong Li). 
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Fig. 3 Ami5 plotted against the relative 
radial distance rsN/-R25- The filled circles 

show 26 SNc la in E/SO galaxies. The open 
circles represent 70 SNe la in spiral galax- 
ies. 



Fig. 4 Distribution of the normalized ra- 
dial distance rsN/-R25 for SNe la in E/SO 
galaxies (solid line) and in spiral galaxies 
(dotted line). 



3.2 Distribution of rsN/^25 

The distribution of SNe la in host galaxies can provide important clues to their progenitors 
and give constraints on the explosion models. In Figure 3, the decline rate Ami5 is plotted 
against the normalized radial distance rsN/-R25- The filled circles mark the 26 SNe la in E/SO 
galaxies, and the open circles, the 70 SNe la in spirals, with the observational error bars shown. 
In the figure SN 2002cs is a slow decliner (with Ami5=1.00) in the elliptical galaxy NGC 6702, 
and SN 2002dj is another slow decliner (with Ami5=1.05) in the elliptical galaxy NGC 5018. 
Moreover, SN 1998es is a SN IQQlT-hke object in a lenticular galaxy. Its host galaxy NGC 632 
has a nuclear starburst with a diameter of 18 kpc, and the position of SN 1998es. 

Figure 4 shows the distribution of the normalized radial distance rsN/^25- Wang (2002) 
first noticed that SNe la might be more centrally concentrated in spiral galaxies than in E/SO 
galaxies. In our study, the fraction of SNe la in E/SO galaxies located in the rsN/R25 > 1-0 
region is about 23%, while this fraction is only 4% for the SNc la in spirals. In other words, 
SNe la in spiral galaxies could be more concentrated towards the optical i?25 radius of their 
host galaxies than in E/SO galaxies. A K-S test shows there is only a 33% probability that 
the ancestors of SNe la in the E/SO and the spiral galaxies have the same radial distribution. 
Furthermore, the plot displays a deficiency of SNc la in the central rsN/-R25 < 0.2 region. The 
data on the SNe la radial distribution closer to center the galaxies can be affected by absorption 
(Windhorst et al. 2002). Note that this effect on the SNe la distribution closer to the center 
of the galaxies may be more serious in spiral galaxies with abundant gas and dust. Therefore, 
SNe la in spiral galaxies may be more centrally located than the above number suggests. 
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3.3 Mni — Ami5 Relation 

The mass of ^^Ni produced during a SN la explosion is the primary physical parameter deter- 
mining the luminosity of the event. One simple method to estimate the synthesized ^^Ni mass 

is a detailed calculation on the assumption that the maximum luminosity of a SN la is propor- 
tional to the instantaneous rate of radioactive decay (Arnett 1982; Arnett et al. 1985: Branch 
1992). Thus the implied bolometric luminosity at light maximum (the luminosity integrated 
from ultraviolet to infrared) can be expressed as 

Lho\ = aS{tn)M^i, (1) 

where Mni is the total mass of ^^Ni produced in the explosion, and a is the ratio of bolometric 

to radioactivity luminosity. S{rYi ) is tlu^ radioactivity luminosity per solar mass of ^^Ni from 
its decay to ^^Co and the subsequent decay to '"'''Fc (e.g., Nadyozhin 1994): 

S = 6.45 X lo43e-*a/8-8 + 1.45 x io43e-*H/iii-3 ^^.g g-i j^-i^ (2) 

where is the rise time in days from the SN explosion to the light maximum. A typical rise 
time of 19 ± 2 days is assumed for SNe la in our analysis (e.g., Branch 1992). 

Prom the analytic solutions for the early-epoch SN la light curves, Arnett (1982) found 
that a should be 1 exactly. On the other hand, according to an early survey of SN la lightcurve 
calculations. Branch (1992) found that a is slightly larger than unity. He concluded a = 1.2±0.2 
and noted that the value of a is independent of the rise time. We assume a = 1 in our analysis 
because the nickel produced above the photosphere may not make significant contribution to the 
luminosity (e.g., Nugent et al. 1995). Stritzinger & Leibundgut (2005) noted that light curves 
of two 3-D SNe la explosion models recently synthesized by the MPA group are consistent with 
a = 1. Using Arnett's rule (a = 1.0) which states that the maximum luminosity of a SN la 
is equal to the instantaneous energy deposition rate from the radioactive decays within the 
expanding ejecta (Arnett 1982; Arnett et al. 1985), and assuming a rise time to bolometric 
maximum of 19 days, one can derive from Equations (1) and (2) a simple formula for the 
maximum luminosity of a SN la (e.g., Stritzinger & Leibundgut 2005): 

L^,^ = (2.0±0.3)xl0«(^)ergs-i. (3) 

The uncertainty of the derived ^^Ni mass is primarily determined by the errors in the adopted 
bolometric rise time and the transformation efficiency of the trapped 7 energy. 

It is generally believed that the ^^Ni mass can be better measured by the bolometric light 
curves constructed from the light curves in the ultraviolet, optical and infrared wavelengths. 
However, only a few SNe la have good observations covering from UV to IR, which could be 
used as a training sample for ^^Ni mass determination of other SNe la. In Figure 5, the ^^Ni 
masses derived from 22 well-observed SNe la with UBVRI and infrared light curves, are plotted 
vs. Amis. For consistency in the Amis determinations, the Amis values of the 22 training SNe 
la are taken from Table 1 in Wang et al. (2006). Of the 22 SNe la with '^'^Ni estimates, 17 arc 
from Table 2 in Stritzinger et al. (2006b) and five (SNe 1995E, 1998aq, 1998de, 1999ac and 
1999dq) from Table 1 in Stritzinger et al. (2006a). To derive the ^^Ni masses Stritzinger et al. 
(2006 a, b) made use of Equation (3). In Figure 5, the solid line corresponding to a least squares 
fit of the data yields the following zero points and dispersion: 



Mm = (-0.61 ±0.09)Amis + (1.30 ±0.11) (7 = 0.13 iV = 22, 



(4) 
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Fig. 5 ^^Ni mass derived from maximum luminosity vs. Ami5 for 22 well-observed 

SNe la. The solid line marks the least squares fit. The dotted lines correspond to a 
vertical deviation of 0.13 (i.e., a quoted dispersion "a") about the solid line. Note 
that the solid line in the panel is similar to the relation obtained by Phillips. 



where the quoted dispersion "a" is the simple root mean square (rms) deviation of the points 
about the fit. Our fitting result is very similar to that of Mazzali et al. (2007). Note that the 
least square fit line is similar to the Phillips Mb-Atoi5 relation. A vertical deviation of 0.13 Mq 
(dotted lines) about the solid line is also shown. It is suggested that the maximum luminosity of 
a SN la (Amis) could be a direct consequence of the '''^Ni mass via modelling the light curves of 
SNe la (Arnett 1982; Cappellaro et al. 1997; Mazzah et al. 2001; Kasen 2006). The correlation 
between the ^^Ni mass and Ami5 shown in Figure 5 may be taken as observational evidence 
for this assumption. 

3.4 Distribution of ^^Ni Masses 

The amount of ^^Ni synthesized from burning to nuclear statistical equilibrium (NSE) in SN la 
explosions is probably the most important physical parameter underlying the observed corre- 
lation of SN la luminosities with their light curves. With the ®^Ni mass one can directly study 
the most sensitive part of the explosion and provide constraints on the explosion mechanisms 
(see the review Leibundgut 2000). Here, according to the empirical relation between the ^^Ni 



8 



B. Wang et al. 



60 r 



40 



Spiral 



20 



03fg 



0.0 0.4 8 



1.2 



20 



- 91bg-like ; 



0.0 



91T-like 

V7, 



0.4 0.8 



03fg 



1.2 



Fig. 6 Distribution of the ^^Ni mass for 
28 SNe la in E/SO galaxies (dotted line) 
and 76 in spiral galaxies (solid line) except 
for SN 2001ay which is truly peculiar. The 
black area represents SN 2003fg. 



Fig. 7 Distribution of the ^^Ni mass for 
different SN la types. The dotted line 
shows the 90 normal SNe la. The gray- 
colored area and the single shaded area 
are for 6 SN 1991bg-like events and 13 SN 
1991T-likc objects, respectively. The black 
area represents SN 2003fg. 



mass and Ami5, i.e., Equation (4), we can estimate the '''^Ni masses and give the distribution 
of the ^^Ni mass for a large number of SNe la (see Figs. 6 and 7). 

As shown in Figures 6 and 7, the ^^Ni mass synthesized in SN la explosions varies from 
0.1 to 1.3M0, with mean error O.17M0. Figure 6 displays the distribution of the ^^Ni mass 
for SNe la in spiral and E/SO galaxies, the black area represents SN 2003fg. To produce SN 
2003fg's unusual luminosity, Howell et al. (2006) calculated, using Equation (1), ^^Ni mass, 
Mni = 1.29 ± O.O7M0, a value that was further confirmed by Jeffery et al. (2006). In Figure 
7 we show the distribution of the ^^Ni mass for different SN la types, and it is indicated that 
different SN la types (i.e., normal SNe la, SN 1991bg-like events and SN 1991T-like objects) 
require the existence of a range in the ^^Ni masses in SN la explosions. 



4 DISCUSSION AND CONCLUSIONS 

As we know above (see Section 3.1), the brighter SNe la prefer to occur in spiral galaxies and 
the fainter ones prefer to occur in E/SO galaxies. However, SNe 2002cs and 2002dj are slow 
decliners in elliptical galaxies. This contradicts with previous conclusion that elliptical galaxies 
produce only fainter SNe la. These two cases probably suggest that the age may not be the 
exclusive factor underlying SN la diversities. Moreover, SN 1998es is a slow decliner in a nuclear 
starburst galaxy, which can probably be taken as an instance of brighter events occurring in 
systems with ongoing star formation (Sullivan et al. 2006). Figure 4 shows that SNe la might be 
more concentrated towards the optical i?25 radius of their host galaxies in spirals than in E/SO 
galaxies. This suggests that SNe la in spiral and E/SO galaxies might have different progenitor 
properties. 



Distribution of ^^Ni Yields 



9 



There is increasing evidence that a single parameter (Amis) cannot account for the diversity 
in SN la (e.g., Bcnctti et al. 2004). As shown in Figure 5, there is a relatively large scatter ~ 
0.13 Mq in the "''^Ni mass- Amis relation. The dispersion may be partly due to uncertainties 
in the distance and reddening adopted in the calculations, but an intrinsic dispersion cannot 
be ruled out. Although it is believed that the ^^Ni mass in SN la explosions is the dominant 
physical parameter affecting both the brightness and light curve shape (or Amis), Woosley 
et al. (2006) proposed that some parameters other than Mni also affect the light curves to 
some extent, such as the total burned mass, the stable iron mass and the degree of direct ^^Ni 
mixing. However, since the empirical MNi-Amis relation is similar to the Phillips Ms-Amis 
relation and the larger sample used in the analysis helps to reduce the effect of random error, 
it is reliable to use the relation to roughly estimate the ^®Ni mass. 

Based on the distribution of the ^^Ni masses above, it is shown that the ^®Ni mass varies 
by at least a factor of ten (e.g., Mni = 0.1 — I.SMq). The explosions of SNe la might produce a 
range in the amount of ^^Ni synthesized from ~0.1 Mq associated with the subluminous objects 
to ^1.3 M0 for the most luminous events. The wide distribution of the '^^Ni masses raises the 
question as to what physical mechanism(s) can account for this dispersion. During the past 
years, a number of potential parameters in the Chandrasekhar-mass deflagration scenario were 
investigated to answer this question, e.g., C-0 ratio, overall mctallicity and central density. 
However, 3D simulations by Ropke & Hillebrandt (2004) indicated that different C-0 ratios 
have a negligible effect on the amount of ^^Ni produced. Also Ropke et al. (2005b) showed 
that there exits a ^ 20% change in the resulting ^^Ni masses via altering the metallicity, 
which is consistent with the analytic prediction of Timmes et al. (2003). Moreover, it has been 
shown that changing the central density of the white dwarf does influence the robustness of the 
explosion. Again, Stritzinger et al. (2006a) suggested that it is unrealistic that any one of these 
parameters, or even a combination of all three, can account for a factor of ten range in the ^^Ni 
mass. 

Besides these parameters, the explosion mechanism itself is more likely to influence the 
^^Ni masses synthesized in SN la explosions (see Stritzinger & Leibundgut 2005). The standard 
Chandrasekhar-mass model is the most popular explosion scenario for SNe la. However, different 
Chandrasekhar-mass models can produce different amount of '"'^Ni which is usually in the range 
of 0.4 - 0.8 (e.g., Mazzali et al. 2001), thus it is an observational fact that a single class 
of Chandrasekhar-mass model cannot account for the explosion properties of all the SNe la. 
Recently, Mazzali ct al. (2007) accomplished a systematic spectral analysis of a large sample 
of well-observed SNe la and mapped the velocity distribution of the main products of nuclear 
burning in SN la explosions. They indicated that a single explosion scenario, possibly a delayed 
detonation (DD) of Chandrasekhar-mass model, may explain most SNe la. From the distribution 
of the ^^Ni mass displayed above, one can find that the ^^Ni mass fall mostly in the range of 
0.4 0.8 M0, thus we suggest that a Chandrasekhar-mass model might account for normal SNe 
la which have the ^^Ni masses from 0.4 Mq to 0.8 Mq. Moreover, Sub-Chandrasekhar WDs 
are plausible candidates for the subluminous objects. We note that one appealing advantage 
offered by sub-Chandrasckhar model is the ability to obtain the progenitor statistics predicted 
by population synthesis calculations (Livio 2000). In addition, the sub-Chandrasekhar model 
can easily explain the observed diversity of SNe la by a one parameter sequence in terms of 
the WD mass (Ruiz-Lapuente et al. 1995), and we also suggested flexible progenitors for SN 
1991bg-nke events with Mni < 0.2 Mq (Cappellaro et al. 1997; Stritzinger et al. 2006a, b). 
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However, it must be noted that the sub-Chandrasekhar model has difficulty in matching the 
observed light curves and spectroscopy (e.g., Hoflich & Khokhlov 1996). 

As shown in Figure?, the ^®Ni mass of SN 2003fg exceeds others remarkably. Hillcbrandt 
et al. (2007) argued that a lop-sided explosion of a Chandrasekhar-mass WD could provide a 
better explanation for SN 2003fg. They suggested that such objects must be rare, since both 
a moderately high ■''^Ni mass (> 0.9 M0) and rather special viewing direction are rare, since 
both a moderately high ^^Ni mass (> 0.9 M©) and a rather special viewing direction are rare. 
However, it is commonly believed that SN 2003fg could be from a super-Chandrasekhar WD 
explosion (e.g., Uenishi et al. 2003; Yoon & Langer 2005; Howell et al. 2006; Jeffery et al. 2006). 
It is proposed that the influence of rotation of the accreting WDs may cause the mass to exceed 
1.4 Mq. Mazzali et al. (2007) also noted that some extremely luminous SNe la may come from 
very rapidly rotating WDs whose mass exceeds Chandrasekhar mass, but these are rare. 

In summary, the wide distribution of the ^^Ni mass suggests that there are probably mul- 
tiple origins of SN la explosion and/or progenitor systems. We argue that sub-Chandrasekhar 
WDs could be the progenitors of some extremely nickel-poor SNe la (with M^i < 0.2 Mq), 
and that super-Chandrasekhar WDs the progenitors of some rare nickel-rich SNe la (with 
Mni > 0.8 Mq). The Chandrasekhar-mass model remains the main scenario for SN la explo- 
sions with a range in ^^Ni masses from 0.4 Mq to 0.8 Mq. To set further constraints on SN 
la explosion/progenitor models, large samples of SNe la with well-observed light curves and 
spectroscopy in nearby galaxies are required to establish the connection of SN la properties 
with the stellar environments of their host galaxies. 
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